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Abstract 

We present global solutions of optically thin, two-temperature black hole accretion disks incorporating 
magnetic fields. We assume that the n7(^-component of the Maxwell stress is proportional to the total 
pressure, and prescribe the radial dependence of the magnetic flux advection rate in order to complete the 
set of basic equations. We obtained the magnetically supported (low-/3) disk solutions, whose luminosity 
exceeds the maximum luminosity for an advection-dominated accretion flow (ADAF), _L^0.4Q;^iEdd, where 
^Edd is the Eddington luminosity. The accretion flow is composed of the outer ADAF, a luminous hot 
accretion flow (LHAF) inside the transition layer from the outer ADAF to the low-/? disk, the low-/3 disk, 
and the inner ADAF. The low-/? disk region becomes wider as the mass accretion rate increases further. In 
the low-/? disk, the magnetic heating balances the radiative cooling, and the electron temperature decreases 
from ^ to ^ lO^K as the luminosity increases. These results are consistent with the anti-correlation 

between the energy cutoff in X-ray spectra (hence the electron temperature) and the luminosity when 
L ^ O.lLEdd, observed in the bright/hard state during the bright hard-to-soft transitions of transient 
outbursts in galactic black hole candidates. 

Key words: accretion, accretion disks — black hole physics — magnetic fields — radiation mechanisms: 
non-thermal 



1. Introduction 

The spectral behavior of black hole candidates (BHCs) 
contains valuable information on physical state of accre- 
tion disks. It is known that galactic BHCs show X-ray 
spectral state transitions in their transient outbursts. The 
system typically undergoes a transition from the low/hard 
state to the high/soft state (so-called hard-to-soft transi- 
tion). In the low/hard state, the X-ray spectrum is ap- 
proximated by a hard power law (photon index ^1.7) with 
an exponential cutoff (E'cut ^ 200 keV) , and the luminos- 
ity is low. In the high/soft state, the X-ray spectrum is 
dominated by a blackbody component emitted from the 
standard accretion disk (Shakura & Sunyaev 1973), and 
the luminosity is high. 

X-ray observations of BHCs have identified two 
types of hard-to-soft transition (e.g., Belloni et al. 
2006; Gierlihski & Newton 2006). One is the bright 
hard-to-soft transition whose transition luminosity is 
high (~ O.S^Edd), where ^Edd = AttcGAI/kcs ~ 1-47 x 
1039 (Af/lOMo) (Kes/0.34 cm2 g-i)"^ erg s"^ is the 
Eddington luminosity, M is the black hole mass, and Kes 
is the electron scattering opacity. The other is the dark 
hard-to-soft transition whose transition luminosity is low 



(< O.lLEdd). 

In the bright hard-to-soft transition, several additional 
X-ray spectral states are identified during the transi- 
tion from the initial low/hard state to the high/soft 
state: bright/hard state (or simply, the brightening of 
the hard state), intermediate state, and very high/steep 
power-law state (e.g., Homan & Belloni 2005; Miyakawa 
et al. 2008). The X-ray spectrum in the bright/hard 
state is approximated by a hard power law as in the 
low/hard state. However, the energy cutoff decreases 
from ^ 200 keV to ^ 50 keV as the luminosity increases 
(in other word, strongly anti-correlates with the lumi- 
nosity) in the bright/hard state while it is roughly con- 
stant (or very weakly anti-correlates with the luminosity) 
around 200 keV in the low/hard state (e.g. Joinet et al. 
2008; Miyakawa et al. 2008; Motta et al. 2009). In gen- 
eral, the energy cutoff is related to the temperature of 
the thermal Comptonizing electrons in an optically thin 
disk (or corona) close to the black hole. Spectral analy- 
ses with thermal Comptonization models (e.g. COMPST 
model in XSPEC introduced by Sunyaev & Titarchuk 
1980, COMPTT model by Titarchuk 1994, and COMPPS 
model by Poutanen & Svensson 1996) also show the anti- 
correlation between the electron temperature and the lu- 
minosity (e.g. Joinet et al. 2008; Miyakawa et al. 2008). 
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In order to account for this anti-correlation, Miyakawa 
et al. (2008) proposed a possible scenario that the heating 
balances the radiative cooling due to the inverse Compton 
scattering; that is, the accretion flow is radiatively efficient 
in the bright/hard state. 

In addition to the appearance of these additional states, 
recent analyses on the X-ray and radio data of BHCs (e.g.. 
Fender et al. 2009) suggest that episodic ejections of rel- 
ativistic jets are associated with the bright hard-to-soft 
transition. The episodic ejections are observed at near 
the peak luminosity of each outburst during the transi- 
tion from the bright/hard state to the very high/steep 
power-law state. (Note that they are not observed in the 
bright/hard state). Therefore, understanding the physical 
mechanism in the bright/hard state may play an impor- 
tant role in understanding the mechanism of the disk-jet 
coupling in BHCs. In this paper, we focus on the physical 
state of the accretion disk in the bright/hard state. 

The standard disk model (Shakura & Sunyaev 1973) has 
been widely and successfully used to account for the black- 
body component. However, optically thin, hot accretion 
disks have been studied to account for hard X-rays from 
BHCs (Thorne & Price 1975; Shibazaki & Hoshi 1975). 
Eardley et al. (1975) and Shapiro et al. (1976) constructed 
a model for an optically thin, two-temperature accretion 
disk in which ions are much hotter than electrons and 
the viscous heating balances the radiative cooling. This 
model, however, is thermally unstable. 

Advection-dominated accretion flows (ADAFs) or ra- 
diatively inefficient accretion flows (RIAFs) were intro- 
duced by Ichimaru (1977) and have been studied exten- 
sively by Narayan and Yi (1994, 1995) and Abramowicz 
et al. (1995). Esin et al. (1997, 1998) found that the max- 
imum luminosity for the ADAF/RIAF is L ^ 0.4Q!^LEdd, 
where a is the viscous parameter (Shakura & Sunyaev 
1973), and showed that the electron temperature is Tc ^ 
lO^-^K. Therefore, the ADAF/RIAF model can account 
for the high-energy cutoff and the low luminosity in the 
low/hard state. However, this model cannot account for 
the relatively low-energy cutoff (i?cut ^ 200 kcV, hence 
the relatively low electron temperature) and the high lu- 
minosity (L ^ O.liEdd) observed in the bright/hard state. 

Luminous hot accretion flows (LHAFs) were proposed 
by Yuan (2001, 2003), in which heat advection works as an 
effective heating and balances the radiative cooling above 
the maximum mass accretion rate for the ADAF/RIAF. 
Although this model is thermally unstable. Yuan et al. 
(2003) concluded that the thermal instability will have no 
effect on the dynamics of the LHAF because the accretion 
timescale is shorter than the growth timescale of the local 
thermal perturbation at high mass accretion rates. The 
LHAF model can partly account for the bright /hard state 
observed in the range of a relatively low luminosity and 
a relatively high energy cutoff (e.g.. Yuan & Zdziarski 
2004; Yuan et al. 2007). However, this model also cannot 
account for the bright /hard state observed in the range of 
a higher luminosity and the lower energy cutoff because 
the electron temperature is a little too high (To ^ lO^K). 
In addition, the anti-correlation between the luminosity 



and the electron temperature is weak. 

In these models, magnetic fields are not considered ex- 
plicitly, and the plasma /? (=Pgas/Pmag) is given as a con- 
stant parameter in general (typically, /? ^ 1). A robust 
mechanism of excitation of magnetic turbulence in accre- 
tion disks is thought to be the magneto-rotational insta- 
bility (MRI), since Balbus and Hawley (1991) pointed out 
its importance. Many local and global magnetohydrody- 
namic (MHD) simulations have investigated the growth 
and the saturation level of the MRI in accretion disks. 
These simulations revealed that the MRI can excite and 
maintain magnetic turbulence and that the Maxwell stress 
generated by the MRI can efficiently transport the angu- 
lar momentum of the disk gas. In addition, several MHD 
simulations and analytical studies suggest that magnetic 
turbulence driven by the MRI can survive even when the 
magnetic pressure is dominant; therefore, highly mag- 
netized accretion disks are astrophysically viable (e.g., 
Shibata et al. 1990; Pessah & Psafiis 2005; Machida et 
al. 2006; Johansen & Levin 2008). 

Machida et al. (2006) demonstrated transitions from 
an ADAF/RIAF- like disk (optically thin, geometrically 
thick, radiatively inefficient, hot, gas pressure dominant 
disk; /3 ~ 5) to a low-/3 disk (optically thin, geometri- 
cally moderately thick, radiatively efficient, cool, mag- 
netic pressure dominant disk; j3 ~ 0.1) by global three- 
dimensional MHD simulations incorporating the radiative 
cooling (see figure 1). When the mass- accretion rate ex- 
ceeds the threshold for the onset of a cooling instability, 
the initial ADAF/RIAF-like disk rapidly shrinks in the 
vertical direction due to the cooling instability. During 
the transition, azimuthal magnetic fluxes inside the disk 
are almost conserved because the timescale of the cool- 
ing instability is shorter than that of the buoyant escape 
of magnetic fluxes from the disk surface. In this way, 
the magnetic pressure becomes dominant and supports 
the disk. Johansen and Levin (2008) performed verti- 
cally stratified shearing box simulations of a local patch of 
such highly magnetized disks. They showed that the MRI 
still survives even in such magnetic pressure dominant 
disks; thus, magnetic fields are still turbulent. Although 
the strong magnetic field reduces the growth rate of the 
Parker instability (Parker 1966) and the MRI, the gener- 
ation of azimuthal magnetic fluxes around the equatorial 
plane still balances the buoyant escape of magnetic fluxes 
from the disk surface. Hence the system can stay in a 
quasi-steady state. We note that such \ow-/3 disks are es- 
sentially different from magnetically-dominated accretion 
flows (MDAFs; Meier 2005) which appear in the inner- 
most plunging region of optically thin disks (see also flg- 
ure 1. We will discuss this issue in the discussion section). 
We focus on the low-/3 disk in this paper. 

Mineshige et al. (1995) suggested that an optically 
thin, magnetic pressure dominated disk emits hard X- 
rays. Pariev et al. (2003) developed a local analytical 
model of an optically thick, geometrically thin, strongly 
magnetized disk which produces spectra quite similar to 
those of the standard disk model. The property of such 
magnetically dominated disks was extensively examined 
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by Begelman and Pringle (2007). Bu et al. (2009) pre- 
sented self-similar solutions of magnetized ADAF/RIAF. 

Oda et al. (2007) constructed an one-temperature 
plasma model of an optically thin accretion disk incor- 
porating magnetic fields on the basis of the results of 
three-dimensional MHD simulations. They assumed that 
the ti7(^-component of the stress tensor is proportional 
to the total pressure, and prescribed the advection rate 
of azimuthal magnetic fluxes in order to complete the 
set of basic equations of the vertically integrated, one- 
dimensional accretion flow in steady state. Oda et al. 
(2009) extended the model to an optically thick disk model 
and Oda et al. (2010) extended it to a two-temperature 
plasma model. They obtained local thermal equilibrium 
solutions, and found a new, thermally stable, low-/? disk 
solution in the optically thin and thick regime. The local 
thermal equilibrium solution of the optically thin low-/3 
disk exists above the maximum mass accretion rate for 
the ADAF/RIAF, and the electron temperature is lower 
than that in the ADAF/RIAF. They concluded that the 
optically thin, low-/3 disk can account for the bright/hard 
state during the bright hard-to-soft transition of BHCs. 
However, the results reported by Oda et al. (2009, 2010) 
were based on the local models in the sense that the inter- 
actions between adjacent annuli of the disk were neglected 
(or, put mathematically, the derivative terms in the basic 
equations were parametrized). Therefore, they could not 
investigate global structures, in particular, composed of 
different types of flows. In addition, accretion flows are 
generally thought to be transonic around the black hole. 
However, the local solution can deviate from the transonic 
solution, in particular, in the inner region where most of 
the X-ray is emitted. 

In this paper, we consider global structures of optically 
thin, two-temperature, black hole accretion disks incor- 
porating magnetic fields. The main purpose is to account 
for the bright/hard state during the bright hard-to-soft 
transition. In particular, we focus on the transition from 
the low/hard state to the bright/hard state which can be 
explained by the transition from the ADAF/RIAF to the 
low-/3 disk. The basic equations are described in section 
2. In section 3 and section 4, we present results of the 
global solutions. Section 5 is devoted to discussion. We 
summarize the paper in section 6. 

2. Model and Assumptions 

2.1. Basic Equations 

In this section, we derive the basic equations for verti- 
cally integrated, one-dimensional steady-state, optically 
thin, two-temperature black hole accretion flows (e.g., 
Kato et al. 2008) incorporating magnetic fields (see also 
Oda et al. 2010) from the resistive MHD equations. We 
adopt cylindrical coordinates {w,ip^z). General relativis- 
tic effects are simulated using the pseudo-Newtonian po- 
tential -0 = -GM/{r - rs) (Paczyhsky & Wiita 1980), 
where G is the gravitational constant, M is the black 
hole mass (we assume M = IOMq in this paper), r = 
{uj^ + z^/^, and rs = 2G'A//c2~2.95xl06(M/10Afo) cm 



is the Schwarzschild radius. For simplicity, the gas is as- 
sumed to consist of protons (ions) and electrons. The 
number density of ions and electrons are equal due to 
charge neutrality, n — Ui = n,,. 
The resistive MHD equations are 

dp 
dt 



V • (pv) = 



(1) 



^ + (..V). 



9(pcec) 
dt 



+ V • [(poEc + Pc) V\-{V- V)po 



ad ' 



dt 



-f V-[(piei+pi)t;]-(t;-V)pi 

= (l-<5hcat)9+-g"= , 



— =VxUxB--w 



(2) 



(3) 



(4) 



(5) 



where p = p\-V p^ is the density, pi = nixn and pc = nion 
are the ion and electron densities, nii and nic are the 
ion and electron masses, v is the velocity, B is the 
magnetic field, j = cV x B/Att is the current density, 
Pgas =Pi+Pe = nk(Ti + Tc) is the gas pressure, pi and 
Pc are the ion and electron gas pressure, Ti and Tg are the 
ion and electron temperature, k is the Boltzmann con- 
stant, ei = (pi/pi) I (7i - 1) and = {p^/ Pe) / he - 1) are 
the internal energy of ions and electrons. Here, 7i = 5/3 
and 7e = 7e (Te) are the specific-heat ratio for ions and 
electrons (details in section 2.4). In the energy equations 
for electrons (3) and ions (4), q+ is the heating rate, q"^ 
is the radiative cooling rate, and is the energy transfer 
rate from ions to electrons via Coulomb collisions. Here, 
^hcat represents the fraction of heating to electrons. In 
this paper, we assume that Jhcat is constant for simplic- 
ity. In the induction equation (5), rjai = c^/Anac is the 
magnetic diffusivity, where CTc is the electric conductivity. 
2.1.1. Azimuthally Averaged Equations 

Three-dimensional global and local MHD simulations 
of black hole accretion disks showed that magnetic fields 
inside the disk are turbulent and dominated by the az- 
imuthal component both in the ADAF/RIAF- like state 
and in the low-/? disk state (e.g., Machida et al. 2006; 
Johansen & Levin 2008). On the basis of results of 
the simulations, we decompose the magnetic fields into 
the mean fields, B = (^0,B^,0), and fluctuating fields, 
SB = {SB^,SB^,6Bz), and also decomposed the velocity 
into the mean velocity, v = (vm^v^^Vz), and the fluctuating 
velocity, 5v = {5vt^,5v^,5vz). We assume that the fluc- 
tuating components vanish when azimuthally averaged, 
l^8v) = {SB) = 0, and that the radial and vertical compo- 
nents of the magnetic fields are negligible compared with 
the azimuthal component, |B^-f-(5B^| 3> |<^Bro|, \SBz\ (see 
the left panel in figure 2). Here, ( ) denotes the azimuthal 
average. 
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We assume that the disk is in a steady state and in hy- 
drostatic balance in the vertical direction. By azimuthally 
averaging equations (1) - (5) and ignoring the second order 
terms of 6v, SB^, and SB^, we obtain 
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(6) 
,(7) 

(8) 
(9) 



where JIko = {GA'I/mY^^ /{m — rs) is the Keplerian an- 
gular velocity. Here, the subscript refers to quantities 
in the equatorial plane. Using the equation of state for 
the ideal gas, the vertically integrated total pressure is 
expressed as 



W„ 



1 



The vertically integrated magnetic tension force is ex- 
pressed as 



-dz 



Airm 



1 2/3-1 
mY+lF^ 
1 2^-1 
ml + p-^ 



Ptotdz 



(17) 



Now we integrate the other basic equations in the ver- 
tical direction. We obtain 



+ 



d V 
am m 
d 

v^T—Pi - VzTrPi = (1 - (Jhcat) - q" 

dm dz 



'Pc)Vz] 

(10) 

'Pi)Vz\ 

(11) 



M - 



-2TTmT,v„ 



,dH-v^) , WtotdlnWu 



dlnm S 91nG7 

Wtot dlnl^Ko 2/3- 



*K0 



m^ 



M(^-40 



-2'Km'^ 



dlnm 14 

(BrnBir_ 



-dz 



0: 



d 



|- [vz{B^)] - i- [v^{B^)] + {V X {5v X 5B))^ 

Oz O'UJ 

-{7?mVx (VxB)}^ , (12) 

where ptot = Pgas + Pmag is the total pressure and Pmag — 
{B'^)/8tt is the azimuthally averaged magnetic pressure. 
The last term on the right-hand side of equation (7) rep- 
resents the magnetic tension force. The third and fourth 
terms on the right-hand side of equation (12) represent 
the dynamo term and the magnetic diffusion term. We 
approximate the induction equation later, based on the 
results of the numerical simulations. 

2.1.2. Vertically Integrated, Azimuthally Averaged 
Equations 

We assume that the radial velocity, v^, the specific 
angular momentum, (. = mv^p, and the plasma /3 (= 
Pgas/Pmag) ai'G independent of z, and that the disks are 
isothermal in the vertical direction for simplicity. The 
surface density, E, the vertically integrated total pressure, 
Wtot, and the half thickness of the disk, H, are defined as 



M 



2'Km'^ mi -|- mo 

ains 



M 



dlnm 



2'Km^ mi 



-ao(Tc 

d\\\H 
d\i\m 

dlnTi 



1 



dlnac(Tc) \ dlnTo 



(18) 

(19) 
(20) 



dlnTc 

'5hcatO+ + Q' 



91nn7 



(21) 



ainS d\nH 



d\nx 



dlnm d In m 

= (l-(5hcat)g+-Q'' 



(22) 



(B^)dz 



[{V X {Sv X SB)}^ 
{^mV X (V X B)}^]drudz + const., 



(23) 



pdz : 



pocxp 



2'KpnH , 
(13) 



Wtnt = 



Ptotdz = 
2TTptotoH , 

Wtot 



Ptoto exp 



1^\ 
2lp) 



(14) 
(15) 



where M is the mass accretion rate (for simplicity, we 
ignore the radial dependence of the mass accretion rate), 
^Ko = ^^^^^Ko is the Keplerian angular momentum and £in 
is the specific angular momentum swallowed by the black 
hole. The second term on the right-hand side of equation 
(19) is a correction resulting from the fact that the radial 
component of the gravitational force changes with height 
(Matsumoto ct al. 1984; Kato et al. 2008). This correction 
is not negligible compared to the pressure gradient force 
in general, and to the effective centrifugal force unless 
H/m <g; 1. We can see this by rewriting this term in the 
form — (Iko'^ /m'^){H/ m)'^{dlnQ}<io/ dlnm) using equation 
(15) . In the energy equations, , ^^'^ 
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the vertically integrated heating rate, radiative cooling 
rate, and energy transfer rate from ions to electrons via 
Coulomb collisions, and ac(Te) = l/[7e(To) — 1] and ai = 
l/(7i — 1). In equation (23), $ is the radial advection rate 
of the azimuthal magnetic flux (hereafter we call it the 
magnetic flux advection rate). 

Wc combine the basic equations as. 



dln{-v^ 



El 
D 



dlnT, 



N2 dlnTi 

15 



N3 

aintu :d ' ^^^^ 



dlnw D ' dlnm 

where D, Ni, N2, and A^3 are the functions of zu, v^^, To, 
and Ti (details in appendix 1). 

We integrate these equations from the outer bound- 
ary using the backward Euler method with the Newton- 
Raphson method (appendix 2). We substitute the Runge- 
Kutta method for the backward Euler method only when 
we fail to solve these equations using the backward Euler 
method. We adjust the parameter t^^ so that flows satisfy 
the regularity condition, D = Ni = N2 = N3 = 0, at the ra- 
dius of the critical point (so-called the shooting method). 
Specifically, we regard solutions as satisfying the regular- 
ity condition when the signs of D, A^i, N2, and change 
around the radius of the critical point and D, Ni, N2, 
and A^3 smoothly increase or decrease with decreasing the 
radius. 

2.2. a -Prescription of the Maxwell Stress Tensor 

Global MHD simulations of radiatively inefficient, ac- 
cretion flows (e.g., Hawley & Krolik 2001; Machida et al. 
2006) showed that the ratio of the azimuthally averaged 
Maxwell stress to the sum of the azimuthally averaged 
gas pressure and magnetic pressure is nearly constant 
{as = ~{Br^B^/4:TT)/{pg^s+Pmag) ~ 0.05 - 0.1), except 
in the innermost plunging region near to the black hole. 
Global and local Radiation-MHD simulations of optically 
thick accretion flows also showed such relations between 
the Maxwell stress and the total pressure (e.g., Hirose et 
al. 2006; Ohsuga et al. 2009). On the basis of the sim- 
ulation results, we assume that the azimuthally averaged 
tniyS-component of the Maxwell stress inside the disk is 
proportional to the total (gas and magnetic) pressure. 



An 



-aptot ■ 



Integrating in the vertical direction, we obtain 



{B^B^p) 

^-dz = -aWt, 



(25) 



(26) 



This is one of the key assumptions in this paper. When 
the magnetic pressure is high, the stress can be high, even 
though the gas pressure is low. We can rewrite this rela- 
tion in terms of the kinematic viscosity, v, as 



where 

JIko 9 In 137/ 



A. 



(27) 
(28) 



Cso = y/Pgaso/po is the sound speed, cao = \/2pmago/po 
is the Alfven speed, and f2 is the angular velocity. We 
can roughly estimate v ~ Vturh x 'turb for turbulent vis- 
cosity, where Vturh and ^turb(~ H) are the characteristic 
velocity and scale length of turbulence, respectively. It 
can be generally expected that Wturb ^ Cso in the gas- 
pressure dominant case while Wturb ^ cao in the mag- 
netic pressure dominant case. Therefore, our formulation 

(Wturb ■ 



V 



CsO 



-CMr) is reasonable. 



2.3. Prescription of the Magnetic Flux Advection Rate 

We complete the set of basic equations by prescrib- 
ing the radial distribution of the magnetic flux advection 
rate on the basis of the result of global three-dimensional 
MHD simulations. Performing the integration in the sec- 
ond term of the induction equation (23), we obtain 



v^{B^)dz = -iJroi?o(^J7)v47riJ 



= [dynamo and diffusion terms] -|- const. 



(29) 



where 
Bo(n7) 



fcTi + fcTc 



1/2 



2t:H 



\ 1/2 



1/2 



(30) 



is the mean azimuthal magnetic field in the equatorial 
plane. 

According to the result of the global three-dimensional 
MHD simulation by Machida et al. (2006), the mag- 
netic flux advection rate at each radius is roughly un- 
changed from before and after the transition from the 
ADAF/RIAF-like disk to the low-^ disk. Following this 
result, we adopt the magnetic flux advection rate as the 
parameter in order to complete the set of the basic equa- 
tions. We then need to prescribe the radial dependence of 
the magnetic flux advection rate. The magnetic flux ad- 
vection rate depends on various mechanisms, such as the 
escape of magnetic fluxes due to the magnetic buoyancy, 
the regeneration of azimuthal magnetic fields by the shear 
motion, the generation of magnetic turbulence through 
the MRI, dissipation of magnetic fields due to the mag- 
netic diffusivity, and magnetic reconnection. If the sum 
of the dynamo term and the magnetic diffusion term is 
zero in the whole region, the magnetic fiux advection rate 
is spatially uniform. The global three-dimensional MHD 
simulation performed by Machida et al. (2006) indicated 
that the magnetic advection rate increases with decreas- 
ing radius, specifically, $ cx ro""'^, in the quasi steady state 
as a result of magnetic dynamo and diffusivity processes. 
In this paper, we parametrize the radial dependence of $ 
by introducing a parameter C as follows: 



V'^out 



-c 



(31) 



where $out is the magnetic fiux advection rate at the outer 
boundary w = Wont and a function of a, M, €out — ^in, 
T'out(=Tc,out + 71_out), and ^out (we illustrated the concept 
of the magnetic flux advection in the right panel in figure 
2). The magnetic flux advection rate is spatially uniform 
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in the case that C = 0, and increases with decreasing radius 
in the case that ^ > 0. In this way, we prescribe the mag- 
netic fiux advcction rate at a radius by setting $out and 

To avoid numerical difficulties concerning the choice of 
£in, we set 4ut-An = 0.5^Ko, Tout =0.3 75Tvir, and /3out = 3 
so that $out has the unique value for given a and M. 

Equation (31) is the second key assumption in this pa- 
per. Prescribing the magnetic flux advection rate enables 
the magnetic pressure to increase when the disk temper- 
ature decreases. In contrast, if we prescribe the plasma 
/3 at each radius instead of the magnetic flux advection 
rate, a decrease in temperature results in a decrease in 
magnetic pressure. This is inconsistent with the results of 
three-dimensional MHD simulations (e.g., Machida et al. 
2006). 

2.4- Energy Equations 

In the conventional theory, the viscous heating was ex- 
pressed as ^ t^^vj [dVL / dzu) ^ where t^i^ is the wLp- 
component of the total stress. Three-dimensional MHD 
simulations of accretion disks indicated that the dissipa- 
tion of turbulent magnetic field energy dominates the to- 
tal dissipative heating rate throughout the disk, and is 
expressed as q+ ^ {B^B^/4:Tr)w {dfl/dzu) (e.g., Hirose et 
al. 2006; Machida et al. 2006; Krolik et al. 2007: Hereafter, 
we refer to it as the magnetic heating rate). We em- 
ploy magnetic heating as the heating mechanism inside 
the disk. Following these simulation results and using the 
a-prescription of the Maxwell stress tensor, equation (25), 
we set the vertically integrated heating rate as follows: 



{B^B^^dhin_ 
(?lnn7 



Air 



dz = —aW tnt^ 



(32) 

We assume that the energy transfer from ions to elec- 
trons occurs via Coulomb collisions and use the energy 
transfer rate given by Stepney and Guilbert (1983) 
and Dermer et al. (1991). We consider bremsstrahlung 
(Svensson 1982; Stepney & Guilbert 1983; Narayan & 
Yi 1995), synchrotron (Pacholczyk 1970; Mahadcvan 
et al. 1996; Esin et al. 1996), and Compton cooling by 
bremsstrahlung and synchrotron photons (Dermer et al. 
1991; Narayan & Yi 1995) as cooling processes [see Oda et 
al. (2010) for details]. The vertically integrated radiative 
cooling rate is expressed as 

Qrad = Obr + Qsy + Qbr,C + ,C ■ (33) 

On the left-hand side of the energy equations. 



1 



1 



7c(re)-i 



3A'3(l/^c)+A'i(l/ge) 
A4(l/0c) 



1 



fli 



1 



(34) 

where Kn are modified Bessel function of the second kind 
of the order n, 0c = kTc/{mcC^) is the dimensionless elec- 
tron temperature. The coefficient ao{Tc) varies from 3/2 



in the case of a non-relativistic electrons to 3 in the case of 
a relativistic electrons (e.g., Chandrasekhar 1939; Esin et 
al. 1997). We assumed that ai = 3/2 in the case of a non- 
relativistic ions because the thermal energy of ions never 
exceeds 10% of the ion rest mass energy in our solutions. 

2.5. Outer Boundary Condition 

We imposed the outer boundary condition at Wo^t = 
lOOOrs, 



9/ad,. 



dm 



0.5£ko('!^ = 



= 



(36) 



(37) 



and 



5/, 



ad,: 



dw 



= , (38) 

where /ad.o and /ad,i are the fraction of the heat advection 
to the heating for electrons and ions (so-called advection 
factors) defined as 



Cad A 



ad.i 



-'heat 



g+ + 



/ad.i — 



Q 



ad.i 



1 



(39) 



(40) 



-(5hoat)Q + 

respectively. Again, to avoid numerical difficulties con- 
cerning a choice of fin , we fixed — Iko {'^ — 3rs ) when we 
compute the outer boundary condition so that the outer 
boundary condition is unique for given a and M . 

Under the boundary condition, we obtain an 
ADAF/RIAF-type solution near the outer boundary for 
a low- mass accretion rate and a low-/3 disk- type bound- 
ary solution for a high-mass accretion rate. In the mid- 
dle range of the mass-accretion rates, we obtain several 
types of boundary solutions, specifically, an ADAF/RIAF- 
type, SLE-type, LHAF-type, low-/3 disk-type, and/or 
intermediate-type solutions. In this case, we choose the 
boundary solution having the highest value of the advec- 
tion factor for ions, that is, the boundary solution closest 
to the ADAF/RIAF-type one. 

3. Global Solutions 

We obtained global solutions of optically thin, two- 
temperature black hole accretion disks by numerically in- 
tegrating the basic equations from the outer boundary. 
The free parameters were M, a, ^hcat, and C. 

We chose a = 0.05 as an example of a pressure-gradient- 
driven flow and a = 0.2 as an example of a viscosity-driven 
flow around the radius of the critical point (Matsumoto 
et al. 1984; Narayan et al. 1997; Nakamura et al. 1997). 
Also, the former corresponds to the case of a low transi- 
tion luminosity from the ADAF/RIAF to the low-/3 disk, 
and the latter corresponds to the case of a high transi- 
tion luminosity, because the maximum luminosity for the 
ADAF/RIAF is L - O.ia^LEdd- 
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We chose C = 0-90 as the fiducial value, so that the 
plasma (3 would be roughly uniform in ADAF/RIAF so- 
lutions (details in section 3.3). When C, has a small value 
(e.g., C = 0)7 we cannot obtain low-/3 disk solutions but 
usual ADAF/RIAF and LHAF solutions. In this case, the 
magnetic pressure becomes negligible compared to the gas 
pressure in the inner region. Hence, the equations reduce 
to that of the conventional model. On the other hand, 
when C, is too large (e.g., Q > 1), the magnetic pressure 
always becomes dominant in the inner region, even at low 
mass-accretion rates. 

The fraction of the magnetic heating (Jhoat is a poorly 
constrained parameter. Yuan et al. (2003) suggested that 
<5hoat ^0.5 be required to fit the spectrum of Sgr A* with 
the ADAF/RIAF model. Sharma et al. (2007) performed 
local sharing box simulations of the nonlinear evolution 
of the MRI in a coUisionless plasma incorporating the 
pressure anisotropy, and showed that ^heat is a function 

of Tc and Ti (approximately, (Shcat = 1 + 'iy/T-jTc )■ 
However, we found no qualitative difference for any value 
of ^hcat in the low-/? disk solutions (see also Oda et al. 
2010). Even in the ADAF/RIAF solutions, we found 
small differences only in the electron temperature and ad- 
vection factors in the inner region. Therefore, we choose 
5\icat = 0.2 as the fiducial value in this paper. 

3.1. Radial Structure and Energy Balance of Global 
Solutions 

First, we show the results for the case a = 0.05, C = 0.9, 
and (5hGat =0.2. Figure 3 shows the radial distribution of 
the electron temperature, Te, the ion temperature, Tj, the 
plasma /3, the magnetic field strength at the equatorial 
plane, Bq, the radial velocity, v^m the surface density, S, 
the ratio of specific angular momentum to Keplerian an- 
gular momentum, ^/^ko? ^-nd the ratio of half thickness 
of the disk to radius, H/tu. Wc denote the radius of the 
critical point at which D ~ Ni— N2 — N^~Qhy the plus 
sign (-I-). The energy balance is illustrated in figure 4. 
The top and second panels show the radial distribution of 
the advection factors for ions /ad.i, and for electrons, /ad,c, 
respectively. In addition, we introduce the total advection 
factor. 



J ad,( 



ad ■ 



Q 



ad.i 



(41) 



in order to classify solutions; this is illustrated in the third 
panel. The bottom panel shows the radial distribution of 
the fraction of the magnetic heating rate to the total heat- 
ing rate for electrons, (5hoat(5'''/(<5hoat'9^ + (5"^), in order to 
illustrate which process mainly heats electrons, the mag- 
netic heating or the energy transfer from ions. The disk 
parameters are summarized in table 1. 

We describe five representative solutions with different 
mass-accretion rates. For an ADAF/RIAF solution (/ad ^ 
0.5 in the entire region) when the mass accretion rate is 
low, A//A/Edd = 1-089 X 10-3 (gray dashed). For a criti- 
cal ADAF/RIAF solution (the minimum value of /ad ~ 0) 
when the mass accretion rate is close to the maximum 



mass accretion rate for the ADAF/RIAF. A//A/Edd — 
8.043 X 10-3 (gray thin sohd). For a LHAF solution (the 
solution having the lowest negative value of /ad, that is, 
the heat advection works as an effective heating most ef- 
ficiently) when the mass accretion rate slightly exceeding 
the maximum mass accretion rate for the ADAF/RIAF, 
Af/AfEdd = 1-224 X 10-2 (long dashed). For a low-/3 disk 
solution (the minimum value of (i ^ Q.l) when the mass 
accretion rate is relatively high, Af/AfEdd = 2.246 x IQ-^ 
(short dashed). For an extremely low-/3 disk solution (the 
minimum value oi j3 0.01) when the mass accretion rate 
is high, Af/AfEdd = 5.984 X 10-^ (solid). 

In the ADAF/RIAF solution, ions are heated by mag- 
netic heating and a substantial fraction of the dissipated 
energy is advected inward. In the outer region, elec- 
trons are mainly heated by energy transfer from ions, 
and heat advection for electrons works as effective coo- 
ing. Meanwhile, electrons in the inner region are mainly 
heated by magnetic heating, and heat advection for elec- 
trons works as effective heating. A substantial fraction 
of the dissipated energy that heats electrons is radiated 
away. 

In the critical ADAF/RIAF solution, the heat advection 
becomes inefficient around GOrg where the energy transfer 
from ions to electrons via Coulomb collisions becomes ef- 
ficient because the surface density increases. As a result, 
the ion temperature slightly decreases, and electrons re- 
ceive almost all of the dissipated energy, which is radiated 
away around this radius. 

In the LHAF solution, the radiative cooling becomes 
efficient in the middle region (20rs ^ w ^ lOOrg) because 
the surface density increases further. Therefore, the elec- 
tron temperature decreases (but slightly). The ion tem- 
perature also decreases because ions are well coupled to 
electrons due to the efficient energy transfer via Coulomb 
collisions. Thus, the heat advection for both electrons and 
ions works as an effective heating in this transition layer 
from the high entropy region to the low entropy region. 
Such heat advection balances the radiative cooling. 

In the I0W-/3 disk solution, the radiative cooling and the 
energy transfer from ions to electrons via Coulomb colli- 
sions become more efficient and overwhelm the heat ad- 
vection except in the innermost plunging region (137 ^ 4rs). 
The gas pressure decreases due to radiative cooling, while 
the magnetic pressure increases due to conservation of the 
magnetic flux advection rate at each radius. As a result, 
the total pressure still remains large. Therefore, the mag- 
netic heating being proportional to the total pressure can 
be large enough to balance the radiative cooling. 

In the extremely low-/3 disk solution, the structure of 
the disk is qualitatively the same as the low-/3 disk solution 
but the temperature is lower and the magnetic pressure 
becomes more dominant. 

We also show the results for the case a = 0.2, Q = 0.9, 
and (Jhcat = 0.2 in figure 5 and figure 6. We obtained 
solutions at mass accretion rates higher than the results 
for the case a = 0.05 basically because the heating rate 
being proportional to a increases. 



8 



H. Oda et al. 



3.2. Force Balance in the Radial Direction 

We describe force balance in the radial direction. The 
non-dimensional pressure gradient force, the effective cen- 
trifugal force, and the magnetic tension force are given by 

Wtot 1 91nm,ot 



E w dlnm 



(42) 



r^2 . 



W M^tot 1 dlnl^Ko 



E w d\nw 



l + (f) 



'\ntu 



and, 



/n 



Wtot 1 2/3-1 



(43) 



(44) 



which are derived from the second term on the left-hand 
side, the first and second terms, and the last term on the 
right-hand side of equation (19), respectively. We note 
that even in the case that i = £ko, the centrifugal force 
has a non-zero (but trivial) value unless H/ ^ 1 due to 
the presence of the correction factor for the gravitational 
force. 

Figure 7 illustrates the radial distribution of the pres- 
sure gradient force (solid), the effective centrifugal force 
(dashed), and the magnetic tension force (dotted) in 
the inner region of the disks (1.6rs < w < lOrg) for five 
representative solutions (the ADAF/RIAF, the critical 
ADAF/RIAF, the LHAF, the low-/? disk, the extremely 
low-/3 disk from top to bottom) for the case a = 0.05, 
C = 0.90, and (5hoat = 0.2. Just outside the radius of the 
critical point, the primary inward force is the effective 
centrifugal force when the mass accretion rate is low. As 
the mass accretion rate increases, the pressure gradient 
force becomes dominant, that is, the flow is pressure- 
gradicnt-driven. The nature of flows very near the radius 
of the critical point is the same as the conventional model 
(Nakamura ct al. 1997; Narayan ct al. 1997). However, 
this is not the case outside this region when the mass 
accretion rate is high; the magnetic tension force is the 
primary inward force. In other word, the pressure gradi- 
ent force pushes the gas inward very near the radius of 
critical point, while the magnetic tension force pushes the 
disk gas inward outside this region. 

Figure 8 illustrates the force balance for the case a = 
0.20, C = 0.90, and (5hcat — 0.2. In contrast to the re- 
sults for the case a = 0.05, the primary inward force is al- 
ways the effective centrifugal force just outside the radius 
of the critical point; that is, the flow is viscosity-driven. 
However, as well as the results for the case a = 0.05, the 
magnetic tension force is the primary inward force outside 
this region when the mass accretion rate is high. 

3.3. Dependence on C 

We investigated the dependence on the parameter C pre- 
scribing the radial dependence of the magnetic flux ad- 
vection rate. We illustrate the radial distribution of the 



plasma /3 and the magnetic flcld strength at the equatorial 
plane, -Bq, in flgure 9, and the electron temperature and 
the ion temperature in flgure 10. The disk parameters are 
a ^ 0.05, Jhoat = 0.2, C = 0.5 (dashed), 0.75 (long dashed), 
0.9 (solid), and 1 (dotted). The three representative so- 
lutions (ADAF/RIAF, LHAF, low-^ disk solutions) are 
illustrated from the top panel to the bottom panel. 

At low mass-accretion rates (i.e., the ADAF/RIAF so- 
lutions), the plasma f3 increases for C < 0.9 and decreases 
for C > 0.9 with decreasing radius, while being roughly 
uniform at ~ 5 for C = 0.9. For this reason, we choose 
C = 0.9 as the flducial value in this paper. 

In the ADAF/RIAF solution, the magnetic heating be- 
ing proportional to the total pressure is insensitive to the 
magnetic pressure because the gas pressure dominates the 
total pressure. In addition, the radiative cooling that is 
contributed by the magnetic field via the synchrotron and 
synchrotron-Compton cooling are inefficient. Therefore, 
the electron and ion temperatures are roughly indepen- 
dent of C, while the plasma /? and Bq strongly depend on 

C- 

As the mass accretion rate increases and exceeds the 
maximum mass accretion rate for the ADAF/RIAF (i.e., 
in the LHAF solutions), the radiative cooling becomes 
efficient in the middle region (20rs ^ w ^ lOOrg). For a 
lower magnetic fiux advection rate (i.e., a smaller value of 
C) , a larger decrease in gas pressure is required in order for 
the magnetic pressure to become high enough to support 
the disk. Therefore, the electron temperature and the ion 
temperature decrease more drastically in the transition 
layer (to ~ 25rs). 

As the mass-accretion rate increases further, the tran- 
sition layer retreats outward and the disk becomes cooler 
and more magnetic pressure dominant. We note that a 
lower magnetic flux advection rate (i.e., a smaller value 
of C) results in lower temperatures and lower plasma /? in 
the inner region. In contrast, the magnetic fleld strength 
attains roughly the same level for different values of C- 

4. Relations between M, L versus Local 
Quantities 

The relations between the mass-accretion rate (or the 
luminosity) and the physical quantities, such as the sur- 
face density, electron temperature at a radius, are widely 
used to understand the X-ray spectral state transition ob- 
served in BHCs. In this section, we consider these rela- 
tions obtained from the global solutions. 

Figure 11 shows the relations between E versus M, Te, 
Ti, and (3 at w = 5rs for the case a = 0.05, Jheat = 0.2, 
C = 0.90 (black), 0.75 (gray), and 0.50 (open diamond). 
We obtained ADAF/RIAF branches in the low mass- 
accretion rate and high-temperature region, and low-/3 
disk branches in the high mass-accretion rate and low- 
temperature region. 

Not the mass-accretion rate, but the luminosity is an 
observable quantity. The electron temperature and the 
Compton y parameter are fundamental parameters of 
X-ray spectral fltting with the thermal Comptonization 



No. ] 

model. To make it easier to compare our results with 
observations, we also show the relation between To and 
the Compton y parameter at w = 5rs versus L in figure 
12. Wc calculated the luminosity L by integrating the 
radiative cooling rate over l.Grg < m < SOOrg for a given 
mass-accretion rate. Note that the most luminous region 
is not the inner boundary, but the region around the ra- 
dius of the critical point where the surface density is high. 
The Compton y parameter is given by 

V = 1 1 + 1 Te, (1 + re,) . (45) 

TTlcC- \ mcC / 

When the luminosity is very low (L ^ lO^^iEdd), both 
the energy transfer from ions to electrons and the radia- 
tive cooling are extremely inefficient, because the surface 
density is very low. Thus heat advection is dominant for 
not only ions, but also electrons, Qad.i^ — ^\ica.t)Q^ and 
Qad,c ~ <5hoatQ^- The clcctrou and ion temperatures are 
determined mainly by these energy equations. Since the 
heat-advection terms and the magnetic-heating term have 
the same dependence on the mass accretion rate, the elec- 
tron and ion temperatures are independent of the mass 
accretion rate (hence the luminosity) , and arc roughly the 
virial temperatures, respectively. 

When the luminosity is below and close to the maxi- 
mum luminosity for the ADAF/RIAF (lO^^LEdd < i < 
0.4a^LEdd ^ O.OOliEdd), the radiative cooling becomes 
efficient for electrons, Qad.o + ^ i^hcatQ"'', because 
the surface density is relatively high. However, the heat 
advection still remains dominant for ions, Qad.i ~ (1 — 
i5heat)<5~''- As a result, the electron temperature weakly 
anti-correlates with the luminosity, while the ion temper- 
ature still remains constant. 

When the luminosity exceeds the maximum luminosity 
for the ADAF/RIAF {L > O.OOlLEdd), the radiative cool- 
ing and the energy transfer from ions to electrons become 
dominant, Q~^^ ^eatQ"*" + and Q'" ^ {1- (5hcat)(3"'"- 
Hence, electrons receive a substantial fraction of the mag- 
netic heating, and the radiative cooling balances the mag- 
netic heating, ~ Q"*". In other words, electrons and 
ions are strongly coupled via Coulomb collisions, and the 
flow is radiatively efficient. As a result, the electron tem- 
perature strongly anti-correlates with the luminosity. 

When the luminosity is high, the Compton y parameter 
correlates with the luminosity for = 0.9 and 0.75, and 
anti-correlates with the luminosity for C = 0.5. 

We show these relations for the case a = 0.2 in figure 
13 and figure 14. In this case, the maximum mass ac- 
cretion rate and the luminosity for the ADAF/RIAF in- 
creases (Mc,A 0.05MEdd and Lc,A ^ 0.016LEdd)- The 
luminosity in the low-/? disk branch exceeds O.l^Edd- 
The Compton y parameter correlates with the luminosity 
for all cases because the electron temperature is relatively 
high and gently decreases with the luminosity compared 
to the results for the case a = 0.05. 
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5. Discussion 

First, we briefly remark why we can obtain low-/3 disk 
solutions. In our model, a decrease in the gas pressure 
results in an increase in the magnetic pressure because due 
to the conservation of the magnetic flux advection rate at 
a certain radius. Therefore, even if an efficient radiative 
cooling decreases the gas pressure, the magnetic pressure 
can increase and support the disk in the vertical direction. 
In addition, the magnetic heating being proportional to 
the total pressure can balance such an efficient radiative 
cooling. In this way, wc can obtain low-/? disk solutions. 

We also remark that such low-/? disks are essentially dif- 
ferent from MDAFs in terms of the energy balance and the 
configuration of magnetic fields. The MDAFs appear in 
the innermost plunging region of optically thin accretion 
disks in global MHD and general relativistic MHD simula- 
tions (e.g.. Fragile & Meier 2009). Outside the innermost 
plunging region, the magnetic fields become turbulent be- 
cause the growth timescale of the MRI is shorter than the 
infiow timescale. The generation of magnetic turbulence 
owing to the MRI balances the dissipation, and the dissi- 
pated energy is converted into thermal energy efficiently. 
In such flows, the magnetic heating rate is consistent with 
the prediction of the a-prescription of the stress tensor, 
that is, proportional to the total pressure. This heating 
balances the heat advection in the ADAF/RIAF and the 
radiative cooling in the low-/3 disk. In addition, the mag- 
netic fields are dominated by the azimuthal component be- 
cause the timescale of the stretching of the magnetic fields 
owing to the shear motion is also shorter than the inflow 
timescale. On the other hand, in the innermost plunging 
region, the ratio of the time scales is reversed, because 
the inflow velocity increases with decreasing the radius 
and exceeds the Alfven velocity. Therefore, the magnetic 
held lines are stretched out in the radial direction before 
the MRI grows, and generates turbulence. Since there is 
no turbulence, no dissipation occurs. As a result, a sub- 
stantial fraction of the gravitational energy is converted 
into the radial infall kinetic energy without being con- 
verted into thermal energy. In such flows, the magnetic 
heating rate predicted by the a-prcscription can no longer 
be valid, and there is almost no heating. Hence, the gas 
pressure and temperature become low, and the flow be- 
comes magnetically-dominated. To sum up, in the low-/? 
disk, the magnetic fields are turbulent and dominated by 
the azimuthal component, and the magnetic heating bal- 
ances the radiative cooling. However, in the MDAF very 
close to the black hole, the magnetic fields are coherent 
and dominated by the radial component, and there is no 
magnetic heating. Although both the \ow-f3 disk and the 
MDAF are cool and magnetically-dominated, they are es- 
sentially different. 

5.1. Transition from ADAF/RIAF to Low-j3 Disk 
through LHAF 

In this subsection, we describe how the ADAF/RIAF 
undergoes a transition to the low-/3 disk as the mass accre- 
tion rate increases. The transition is illustrated schemat- 
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ically in figure 15. 

When the mass-accretion rate is below the maximum 
mass-accretion rate for the ADAF/RIAF, we obtained the 
ADAF/RIAF solutions that are essentially the same as 
the solutions of the conventional model (bottom panel in 
figure 15). 

At the mass-accretion rate slightly exceeding the max- 
imum mass accretion rate for the ADAF/RIAF, the ra- 
diative cooling overwhelms the heat advection, working 
as effective cooling in the middle region {w ^ 50rs). The 
gas pressure decreases due to the radiative cooling while 
the magnetic pressure increases because of the conserva- 
tion of magnetic fiux advection rates given at a radius. 
In this way, the magnetic pressure becomes dominant and 
supports the disk in this region, and the flow undergoes a 
transition from the outer ADAF/RIAF to the low- (3 disk. 
The LHAF appears in this narrow transition layer from 
the outer ADAF/RIAF to the low-/3 disk because such a 
flow configuration results in the negative entropy gradi- 
ent. In the inner region {-nj ^ lOrg), the radiative cool- 
ing is still inefficient because the surface density decreases 
and the radial velocity steeply increases with decreasing 
the radius (such a feature is more prominent for lower 
a; see also Nakamura et al. 1997; Narayan et al. 1997). 
Therefore, the flow returns to the ADAF/RIAF. As a re- 
sult, the flow is composed of the outer ADAF/RIAF, the 
LHAF inside the narrow transition layer, the low-/3 disk, 
and the inner ADAF/RIAF (third panel in figure 15). 

We note that such an inner ADAF/RIAF cannot be ob- 
tained from self-similar solutions of optically thin disks. 
In the self-similar solutions, the radiative cooling can be 
efficient, even in the inner region because the surface den- 
sity increases and the radial velocity gently increases with 
decreasing radius (S cx n7~^/^, oc 137"-'^/^). However, 
the self-similar solutions are no longer valid in the inner 
region because of the transonic nature of the flow. 

As the mass-accretion rate increases further, the radia- 
tive cooling becomes more efficient over the whole region. 
Thus, the transition layer between the outer ADAF/RIAF 
and the low-/3 disk retreats, and the inner ADAF/RIAF 
region diminishes. In other words, the low-/3 region be- 
comes wider. Eventually, the whole region becomes the 
low-/3 disk, except in the innermost ADAF/RIAF region 
around the radius of the critical point (first panel in figure 
15). 

5.2. Dynamical Property 

Transonic flows around black holes with the a- 
prescription of the stress tensor are divided into two 
classes according to the value of a, a pressure-gradient- 
driven flow for small a and a viscous-driven flow for 
large a (e.g., Matsumoto et al. 1984; Nakamura et al. 
1997; Narayan et al. 1997). We calculated the solution 
with a = 0.05 as an example of the pressure-gradient- 
driven flow and the solution with a = 0.2 as an example 
of the viscous-driven flow. 

First, we discuss the case for a ~ 0.05. When the mass- 
accretion rate is around or above the maximum mass- 
accretion rate for the ADAF/RIAF, the surface density 



decreases and the radial velocity increases sharply with 
decreasing the radius in the inner region {w ^ lOrg). Just 
outside the radius of the critical point, the angular mo- 
mentum of the flow approaches and exceeds the Kcplcrian 
angular momentum. In this region, the primary inward 
force is the pressure gradient force. 

Next, we discuss the case for a = 0.2. In this case, the 
radial distribution of the surface density is relatively flat 
compared to the case for a = 0.05 and the angular momen- 
tum is always below the Keplerian angular momentum. 
Just outside the radius of the critical point, the primal 
inward force is the effective centrifugal force; that is, the 
gas falls inward due to losses of its angular momentum. 

Our solutions have essentially the same nature of tran- 
sonic flows just outside the radius of the critical point as 
the conventional model. In addition to this nature, we 
found that the primal inward force is the magnetic ten- 
sion force at some distance from the radius of the critical 
point; that is, the magnetic tension pushes the disk gas 
more strongly than the other forces in this region. This is 
a new finding of our results. The magnetic field can con- 
tribute to not only the vertical hydrostatic balance, but 
also the radial force balance in some part of the disk. 

5.3. Application to Bright/Hard State during Bright 
Hard-to-Soft Transition 

The relations between E versus M, To, Ti, and /3 are 
consistent with the thermal equilibrium solutions pre- 
sented by Oda et al. (2010). We note that the LHAF 
branches do not appear at such an inner region, because 
the transition layer from the outer ADAF/RIAF to the 
low-/? disk appears in the middle region {zu ^ 50rs) first, 
and retreats as the mass accretion rate increases. 

We calculated the luminosity by integrating the radia- 
tive cooing rate over 1.67's < ti? < SOOrg. When the lu- 
minosity is very low (i ^ lO^^iEdd for a = 0.05 and 
L ^ 10~*LEdd for a = 0.2), the electron temperature is 
roughly independent of the luminosity and ~ lO^^'^K. 
This indicates that the cutoff energy in the X-ray spec- 
trum is independent of the luminosity (the clear cutoff at 
such a low luminosity, however, may not be detectable). 

When the luminosity is below the maximum luminos- 
ity for the ADAF/RIAF, {L < O.OOlLsdd for a = 0.05 
and L ^ 0.016iEdd for a = 0.2), the electron temperature 
weakly anti-correlates with the luminosity in the range 
from ^ 10^'^K to ^ 10^^'^K. Note that the minimum elec- 
tron temperature for the ADAF/RIAF is roughly inde- 
pendent of the value of a. This feature agrees with the 
result presented by Esin et al. (1998). This weak anti- 
correlation between the electron temperature and the lu- 
minosity in the high electron temperature region can be 
consistent with the weak anti-correlation between the en- 
ergy cutoff and the luminosity observed in the low/hard 
state. However, the ADAF/RIAF cannot account for the 
electron temperature lower than ~ 10^'^K. 

When the luminosity is above the maximum luminosity 
for the ADAF/RIAF, the electron temperature strongly 
anti-correlates with the luminosity in the range from ~ 



No. ] 

lO^K to ^ This strong anti-correlation in the rela- 

tively low electron temperature and high luminosity region 
can be consistent with the anti-correlation between the en- 
ergy cutoff and the luminosity observed in the bright/hard 
state. Therefore, we conclude that the low-/3 disk can ac- 
count for the bright/hard state, and that the transition 
from the ADAF/RIAF to the low-/3 disk corresponds to 
the transition from the low/hard state to the bright/hard 
state during the bright hard-to-soft transition. 

We also touch on the possibility of the dark hard-to-soft 
transition during which the system immediately under- 
goes a transition from the low/hard state to the high/soft 
state at a low luminosity. When ( has a small value, and 
thus, the magnetic flux advection rate is very low in the 
inner region, we could not obtain an optically thin global 
solution at a mass-accretion rate higher than the maxi- 
mum mass accretion rate for the ADAF/RIAF. In this 
case, we expect that the ADAF/RIAF undergoes a tran- 
sition to an optically thick disk (e.g., the standard disk, 
the slim disk, and an optically thick low-/? disk) with the 
LHAF transition layer. This might correspond to the dark 
hard-to-soft transition. 

5.4- What Mechanism Determines the Magnetic Flux 
Advection Rate? 

Wc propose a possible scenario that what mechanism 
determines whether the system undergoes the bright 
hard-to-soft transition or the dark hard-to-soft transi- 
tion. Three-dimensional MHD/Radiation-MHD simula- 
tions showed that the polarity of the azimuthal magnetic 
field inside the disk can change in time alternately in 
the ADAF/RIAF and the standard disk, that is, the gas 
pressure dominant disk (e.g., Nishikori ct al. 2006; Shi 
et al. 2010). Large filament-like structures of magnetic 
fields emerge from the mid plane of the disk and rise up, 
roughly, in the growth timescale of the Parker instabil- 
ity. Subsequently, reversals of azimuthal magnetic fields 
take place in the mid plane, and the polarity of azimuthal 
magnetic fields alternates successively. (Note that such 
a change in polarity can be suppressed in the low-/3 disk 
because the growth timescale of the Parker instability is 
quite long). The pattern of the polarity is symmetrical 
with respect to the equatorial plane roughly, but not ex- 
actly. Hence, the magnetic fiux advection rate inside the 
disk can change in time. 

If the mass-accretion rate exceeds the threshold for the 
onset of the cooling instability when the magnetic fiux 
advection rate is high on average, the ADAF/RIAF will 
evolve toward the low-/3 disk with the LHAF transition 
layer. In this case, we can expect that the system under- 
goes the bright hard-to-soft transition. On the other hand, 
if the mass accretion rate exceeds the threshold when the 
magnetic flux advection rate is low, the ADAF/RIAF will 
evolve toward an optically thick disk with the LHAF tran- 
sition layer. In this case, we can expect the dark hard-to- 
soft transition. 
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5.5. Beyond the Low-ji Disk: Episodic Ejections of 
Relativistic Jets 

According to Pessah and Psaltis (2005), the MRI is sta- 
bilized for toroidal Alfven speeds exceeding the geometri- 
cal mean of the sound speed and the rotational speed of 
the disk gas (wa ^ \/c^vkq)- In mildly low-/? disk solutions 
at moderately high mass accretion rates, this condition is 
not satisfied, that is, the MRI is not stabilized. On the 
other hand; in extremely low-/? disk solutions at very high- 
mass accretion rates, this condition is satisfied in a certain 
region, and thus, the MRI can be stabilized. Hence, we ex- 
pect that no magnetic turbulence can be driven, and then 
no magnetic heating can occur in this region (although 
the amplification of azimuthal magnetic fields due to the 
shear motion may still survive) . Therefore, the cooling in- 
stability will occur and the disk will shrink further in the 
vertical direction. As a result, the magnetic field strength 
and energy inside the disk can be amplified further. Such 
a drastic and further increase in the magnetic field energy 
may lead to an explosive energy release (e.g, Shibata et 
al. 1990; Yuan et al. 2009). The timing analyses on the 
X-ray spectrum and the radio emission of BHCs suggest 
that episodic ejections of relativistic jets take place dur- 
ing the transition from the bright/hard state to the soft 
state (e.g.. Fender et al. 2009). This explosive energy re- 
lease from the low-/? disk may be an origin of the episodic 
ejections of relativistic jets. 

6. Summary 

We have calculated vertically integrated, one- 
dimensional, steady-state global solutions of optically 
thin, two-temperature, black hole accretion disks in- 
corporating the mean azimuthal magnetic fields. We 
have obtained the magnetic pressure dominant (low-/?) 
disk solutions at high mass-accretion rates. We have 
concluded that the low-/3 disk can account for the 
bright /hard state (or simply, the brightening of the hard 
state) observed during the bright hard-to-soft transition 
in transient outbursts of BHCs. 

We have assumed that the ti7(y9-component of the az- 
imuthally averaged Maxwell stress tensor is proportional 
to the sum of the gas and magnetic pressure. We have also 
prescribed the radial distribution of the magnetic fiux ad- 
vection rate by introducing the parameter Q in order to 
complete the set of basic equations. Accordingly, a de- 
crease in temperature results in an increase in magnetic 
pressure under conservation of the magnetic fiux advec- 
tion rate at each radius, and the magnetic heating being 
proportional to the total pressure can balance the radia- 
tive cooling. 

When the mass-accretion rate is below the maximum 
mass accretion rate for the ADAF/RIAF, we obtained the 
usual ADAF/RIAF solutions. When the mass accretion 
rate is just beyond the maximum mass accretion rate for 
the ADAF/RIAF, we obtained the solutions of the fiow 
composed of the outer ADAF/RIAF, the LHAF inside the 
transition layer from the outer ADAF/RIAF to the low- 



Magnetically Supported Accretion Disk 



12 



H. Oda et al. 



/3 disk, the low-;3 disk, and the inner ADAF/RIAF. This 
low-/3 disk region becomes wider as the mass-accretion 
rate increases further. Eventually, the whole region be- 
comes the low-/3 disk, except in the innermost plunging 
region around the radius of the critical point. The elec- 
tron temperature decreases from ~ 10^'^K to ~ lO^K with 
increasing the luminosity above the maximum luminosity 
for the ADAF/RIAF {L > OAa'^LEdd)- This is consistent 
with the anti-correlation between the energy cutoff in the 
X-ray spectrum (hence the electron temperature) and the 
luminosity when L^O.lLEdd, observed in the bright/hard 
state during the bright hard-to-soft transition of BHCs. 

When we assumed very low magnetic flux advection 
rates, we could not obtain the low-/3 disk solutions. In 
this case, we expect that the flow will be composed of the 
outer ADAF/RIAF, the LHAF inside the transition layer 
from the outer ADAF/RIAF to an inner optically thick 
disk, and the inner optically thick disk at the maximum 
mass accretion rate for the ADAF/RIAF. This might cor- 
respond to the dark hard-to-soft transition. 

Although we could obtain the extremely low-/? disk so- 
lution at a very high mass accretion rate, the MRI can be 
stabilized in such an extremely low-/? plasma. Hence, the 
magnetic heating due to turbulent magnetic fields may 
not occur under such an extremely low-/? regime. In this 
case, the cooling instability will occur and the disk will 
shrink further in the vertical direction. As a result, the 
magnetic field energy can be amplified further. Such a 
drastic increase in magnetic field energy may lead to an 
explosive energy release. This can be an origin of the 
episodic ejections of relativistic jets observed during the 
bright hard-to-soft transition. 
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Appendix 1. Combined Form of Basic Equations 

Combining the basic equations, we rewrite the set of 
the basic equations in the following form: 



Ax — b 



(Al) 



Ai/353-l-# -Air„,9o3 ai-Airi^03 



(A2) 



dlnzu 
91nTe 



dlnzu 
dlnTi 



[Vol. , 



(A3) 



dlnzu 

"KO q '"°"kO nf-| I 1-3"^ 

B(l + p-l) din™ ^^>+ x + 

%t >+A„/3o3/3-^ (^^^y^-Hl-2c) 
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To 
A, 



To 

Ti + Te ' 



^+A.fe/3-i(;iiiI^511 + l-2C) + ;iiiI^J 



(A4) 



fcTi + fcTe Ti 
J i i — 



(rTij + rrio) 



_ 2 + 5/?-i 

— n , 00-1 ' P03 



2 + 3^- 



2 + 



<^hcat O? 1 



Fo M 2 



1 - ^heat _ 1 

Fi ^i 2 



271^2 g'° - Q-^ 



Sq 

Si 

A,{T,)= a,{T,) ( 1 + 



ad 



M kTo/ (TOi-t-mo) 

M kT-J (mi -I- TOe) 
dlnac(Tc) 
dluTc 

We define the denominator and the numerators in equa- 
tion (24) as 



L» = detA 



N3 



= adjA • b , 



(A5) 



where adjA is the adjugate matrix of A. 

Appendix 2. Integration Method 

Introducing X = [lni;ro,lnTe,lnri], R = Inzu, and F = 
[Ni/D,N2/D,N3/D], we rewrite equation (24) in the fol- 
lowing form: 



(A6) 



The difference equation in backward Euler method is given 

by 

Xi+i — Xi 



— F{Ri+i,Xi+i) 



(A7) 



where ARi = Ri-^i — Ri. Here we define the residual 4> as 
^{X^+^) = -X^- AR,F{R,+i,X,+i) . (A8) 
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We solve = using Newton-Raphson method. Let Xf_^-^ 
be the current approximation. Then the next approxima- 
tion is given by 

X^+i=Xr+i-[0'(Xr+i)]-'0(Xr+i) . (A9) 
Here (p' denotes the derivative (p with respect to Xi-^-i. 
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ADAF/RIAF ADAF/RIAF 




Fig. 1. Transition from an ADAF/RIAF-like disk to a low-/3 disk. The left panel illustrates the ADAF/RIAF-likc disk state at a 
low mass accretion rate before the transition. The right panel illustrates the low-/3 disk state at a moderately high mass accretion 
rate after the transition. Solid curves depict magnetic filed lines. In both states, the magnetic fields inside the accretion disk are 
turbulent and dominated by the azimuthal component except in the innermost plunging regions (the MDAF regions). 




Fig. 2. Schematic pictures of magnetic field lines inside the accretion disk. Left; Structures of turbulent magnetic fields driven by 
the MRI. Right: Magnetic fiux advcction with mass accretion. 
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Fig. 3. Radial structures of optically tiiin, two-tcmpcraturc accretion disks for the case 

M/MEdd = 2.089 X 10"^ (gray dashed), 

8.043 X 10-3 (gray thin solid), 1.224 X 10~^ (long dashed), 2.246 X 10"^ (short dashed), and 5.984 X 10'^ (solid), which correspond 
to the ADAF/RIAF, the critical ADAF/RIAF, the LHAF, the Iow-/3 disk, and the extremely low- /J disk solutions, respectively. The 
disk parameters are a = 0.05, f = 0.90, and Si^cat = 0.2. The radius of the critical point is denoted by the plus sign (+). 
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Fig. 4. Energy balance for the case a = 0.05, f = 0.90, and <5ij(,at = 0-2. Top: the advcction factor for ions. Second: the advection 
factor for electrons. Third: the total advection factor. Bottom: the fraction of the magnetic heating to the total heating rate for 
electrons. 
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Fig. 5. Radial structures of optically tiiin, two-tcmpcraturc accretion disks for the case 

Af/A>/Edd = 3.993 X 10~^ (gray dashed), 

5.712 X 10~^ (gray thin solid), 1.202 X 10^^ (long dashed), 1.445 X 10~^ (short dashed), and 3.631 X 10"^ (solid), which correspond 
to the ADAF/RIAF, the critical ADAF/RIAF, the LHAF, the Iow-/3 disk, and the extremely Iow-/3 disk solutions, respectively. The 
disk parameters are a = 0.20, f = 0.90, and <5ij(,at = 0-2. 
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(solid), effective centrifugal force, /c, (dashed), and magnetic 
tension force /m, (dotted). 
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Fig. 8. Force balance in the radial direction for the case 
a = 0.20, C = 0.90, and 5heat = 0.2. From top to bot- 
tom. M/MEdd = 3.993 X 10-3, 5.712 x lO^^, 1,202 x IQ-^, 
1.445 X 10-^, and 3.631 X 10-^. Pressure gradient force, /p, 
(solid), effective centrifugal force, /c, (dashed), and magnetic 
tension force /m, (dotted). 
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Fig. 9. Radial distribution of /3 (left) and Bo (right). Three representative solutions (ADAF/RIAF, LHAF, and low-/3 disk) are 
illustrated from top to bottom. The disk parameters are a = 0.05, 5hcat = 0.2, ^ = 0.50 (dashed), 0.75 (long dashed), 0.90 (solid), 
and 1.00 (dotted). 
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Fig. 10. Radial distribution of tlio electron temperature (left) and the ion temperature (right). Three representative solutions 
(ADAF/RIAF, LHAF, and low-/3 disk) are illustrated from top to bottom. The disk parameters axe a = 0.05, 5heat = 0.2, f = 0.50 
(dashed), 0.75 (long dashed), 0.90 (solid), and 1.00 (dotted). 
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Fig. 11. Relation between E versus M (top), Tc (second), Ti (third), and /3 (bottom) at = 5rs. The disk parameters are a = 0.05, 
^hoat = 0.2, C = 0.90 (black), C = 0.75 (gray), and C = 0.50 (open diamond). 
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Fig. 12. Relation between Tg (left), and y (right) at ro = Sr^ versus L. The disk parameters are a = 0.05, ^h^at = 0-2, C = 0.90 
(black), C, = 0.75 (gray), and C, = 0.50 (open diamond). 
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Fig. 13. Relation between S versus M (top) Tc (second), Tj (tliird), and /3 (bottom) at = 5rs. Tiie disk parameters are a = 0.2, 
<5hoat = 0.2, C = 0.90 (black), C = 0.75 (gray), and C = 0.50 (open diamond). 
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Fig. 14. Relation between Tg (left), and y (right) at ro = Sr^ versus L. The disk parameters are a = 0.20, ^h^at = 0-2, C = 0.90 
(black), = 0.75 (gray), and = 0.50 (open diamond). 
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Fig. 15. Schematic pictures of configuration of accretion disks for various mass accretion rates. ADAF/RIAF regions are denoted 
by gray, LHAF regions by white, and low-/3 disk regions by black. The bottom panel shows the low/hard state at a low mass accretion 
rate. The top and second panels show the bright/hard state at high mass accretion rates. The middle panels show intermediate 
states during the transition from the ADAF/RIAF to the low-/3 disk, in which the LHAF appears. 
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Table 1. Disk parameters, luminosity, and radius of critical point 



a 


C 


(^hoat 




M/AfEdd 




L/LEdd 




fcrit/^s 


Type 


0.05 


0.90 


0.2 


1.4083924375 


2.089 X 10' 


-y 


1.430 X 10- 


-4 


2.4443618 


ADAF/RIAF 








1.5674186805 


8.043 X 10- 


-3 


1.450 X 10- 


-3 


2.2445909 


Critical ADAF/RIAF 








1.7231616728 


1.224 X 10- 


-2 


5.079 X 10- 


-3 


2.1117202 


LHAF 








1.7794718322 


2.246 X 10' 


-2 


1.241 X 10- 


'2 


2.0844964 


Low-/3 disk 








1.8124371548 


5.984 X 10- 


-2 


3.856 X 10- 


-2 


2.0794780 


Extremely Low-/? disk 


0.2 


0.90 


0.2 


0.5974897133 


3.993 X 10" 


-3 


2.965 X 10- 


'4 


4.8020671 


ADAF/RIAF 








0.8908277177 


5.712 X 10- 


-2 


1.149 X 10- 


-2 


3.1304441 


Critical ADAF/RIAF 








1.4083890398 


1.202 X 10- 


-1 


9.244 X 10- 


'2 


2.3096968 


LHAF 








1.4445944039 


1.445 X 10- 


-1 


1.209 X 10- 


-I 


2.3014862 


Low-/? disk 








1.5259917051 


3.631 X 10- 


-1 


3.622 X 10- 


-1 


2.3316265 


Extremely Low-/? disk 


0.05 


1.00 


0.2 


1.3491502072 


2.089 X 10- 




1.522 X 10- 


-4 


2.2885442 


ADAF/RIAF 








1.6082077453 


1.546 X 10- 




5.102 X 10- 


-3 


1.9937257 


LHAF 








1.6680056446 


2.353 X 10- 


-2 


1.051 X 10- 


-2 


1.9548674 


Low-/? disk 


0.05 


0.75 


0.2 


1.4389894702 


2.089 X 10- 


-A 


1.298 X 10- 


-4 


2.5295339 


ADAF/RIAF 








1.8135734019 


8.830 X 10- 


-3 


4.700 X 10- 


-3 


2.2610183 


LHAF 








1.8371197818 


2.046 X 10- 


'2 


1.259 X 10- 


-2 


2.2653863 


Low-/? disk 


0.05 


0.50 


0.2 


1.4722433799 


3.162 X 10- 


-3 


1.954 X 10- 


-4 


2.5178226 


ADAF/RIAF 








1.8419696927 


6.992 X 10" 


-3 


4.181 X 10- 


-4 


2.5116158 


LHAF 








1.8441905224 


1.564 X 10- 


-2 


9.605 X 10- 


-3 


2.5183405 


Low-/? disk 



